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The  regularly  sampled,  sleep-wake  behavior  of  a  mouse  can 
be  described  stochastically  as  a  time  series  in  which  a  binary 
valued  random  variable  fluctuates  between  the  two  possible 
states.   It  is  readily  appreciated  that  such  a  series  is  not 
entirely  random,  but  instead  exhibits  significant  local  serial 
correlation  (LSC)  across  a  range  of  time  intervals  (lags). 
This  correlation  results  in  part  from  the  organization  of  the 
alternating  episodes  of  the  two  states  into  consolidated  "bouts" 
of  sleep  and  wake  (1) .   Serial  correlation  beyond  the  mean 
length  of  the  consolidated  bouts,  however,  would  be  of  signifi- 
cant interest  and  could  arise  in  one  of  two  ways;  i)  as  the 
result  of  a  rhythmic  process  governing  the  regular  recurrence 
of  bouts,  or  ii)  as  the  result  of  an  arhythmic  process  in  which 
the  lengths  of  successive  bouts  are  serially  correlated  with 
one  another  and/or  have  skewed  distributions.   Theoretically, 
the  two  processes  should  be  distinguishable  on  the  basis  of 
their  resultant  serial  correlation  functions.   Practically, 
recognition  of  the  differences  in  correlation  requires  time 
series  of  sufficient  length  and  resolution. 

The  majority  of  work  examining  temporal  distribution  in 
mammalian  behavioral  state  data  has  embraced  the  hypothesis 
that  a  short-period  "ultradian"  rhythm  was  the  source  of  ob- 
served LSC.   In  their  original  characterization  of  the  REM 
and  NREM  sleep  states  in  human  subjects,  Dement  and  Kleitman 
observed  a  REM/NREM  alternation  with  a  period  of  approximately 
90  minutes  (2).   The  regularity  of  the  REM/NREM  alternation 
and  observations  of  apparent  local  periodicity  in  other 


psychological  and  physiological  variables  led  Kleitman  to 
postulate  the  existence  of  a  unifying  ultradian  oscillator,  a 
"Basic  Rest-Activity  Cycle"  (BRAC) ,  which  transcended  the 
sleep/wake  boundary  and  influenced  the  expression  of  a  variety 
of  behaviors  in  both  states  (3).   This  attractive  hypothesis 
led  to  a  number  of  attempts  to  demonstrate  ultradian  rhythms 
in  physiology.   Most  of  these  attempts,  however,  have  been 
complicated  by  difficulties  in  (i)  defining  variables  that  can 
be  realiably  measured  across  the  entire  circadian  cycle  and 
that  demonstrate  clear  ultradian  rhythmicity,  (ii)  obtaining 
data  series  of  sufficient  length  and  resolution  to  allow  defin- 
itive statement  about  the  presence  or  absence  of  periodicity, 
and  (iii)  allowing  for  serial  correlation  arising  from  other 
sources  in  the  analysis  of  time  series  data  and,  especially, 
in  the  interpretation  of  spectral  and  autocorrelation  results  (4) 

In  this  paper,  we  report  data  obtained  from  study  of  the 
temporal  sleep-wake  distribution  of  the  mouse.   Three  factors 
make  our  experimental  system  particularly  well  suited  to  study 
of  the  temporal  organization  of  sleep-wake  state  and  the  possible 
contribution  of  a  BRAC.   First,  we  are  currently  capable  of 
collecting  continuous  measurements  of  sleep-wake  state  sufficient 
in  length  (up  to  280  days)  and  resolution  (sampled  every  10 
seconds)  to  allow  rigorous  determination  of  the  significance 
and  consistency  of  apparent  ultradian  and  circadian  rhythmicity. 
Second,  the  alternation  of  sleep  and  wakefulness  has  obvious 
intuitive  validity  as  a  physiological  marker  of  the  rodent  BRAC, 


should  one  exist.   Finally,  though  the  marginal  (first-order) 
probabilities  of  sleep  and  wakefulness  in  the  mouse  demonstrate 
clear  circadian  periodicity,  both  states  are  distributed 
throughout  the  entire  circadian  cycle;  i.e.  the  probability 
P(cf>)   of  sleep  varies  with  phase,  <$> ,  but  almost  never  becomes 
one  or  zero.   Thus,  if  a  BRAC  continues  throughout  the  circadian 
cycle,  its  effect  on  the  expression  of  sleep  and  wakefulness 
and,  in  particular,  on  the  correlation  structure  of  those 
variables  should  be  measurable  at  all  phases  of  the  mouse's 
circadian  cycle.   This  fortuitous  feature  of  mouse  sleep  data 
obviates  the  use  of  questionable  pairings  of  parameters  from 
wakefulness  and  sleep  on  the  assumption  that  they  represent 
multiple  indicators  of  the  same  rhythmic  process  (e.g.  REM 
sleep  and  daydreaming  fantasy  (5) ) . 

Nine  male  mice  (Mus  musculus)  were  studied  while  individually 
housed  in  light-tight,  insulated  enclosures  illuminated  by 
externally  controlled,  water- jacketed  fluorescent  lamps.   For 
this  study,  data  were  collected  under  two  lighting  conditions; 
either  continuous  darkness  ("DD"),  or  12  hours  of  light  followed 
by  12  hours  of  darkness  ( "LD12 : 12" ) .   Food  and  water  were 
available  ad  libitum.   EEG  recording  techniques  and  automated 
state-scoring  methods  were  performed  as  described  elsewhere  (6). 
State  was  determined  every  10  seconds  as  either  active  wake 
(wheel  running),  quiet  wake,  REM  sleep  or  NREM  sleep.   (For  the 
study  reported  here,  active  and  quiet  wake  were  pooled  to 
form  "total  wake",  while  REM  and  NREM  sleep  were  pooled  as 
"total  sleep.")   We  selected  fourteen  sections  from  the  full-length 
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data  trains  for  inclusion  in  the  analysis  on  the  basis  of  (i) 
minimal  interruption  and  lost  data,  (ii)  relative  stationarity 
of  dominant  circadian  period  (for  DD  samples),  and  (iii)  total 
length  between  20  and  30  complete  circadian  cycles.   Eight 
samples  were  selected  from  the  DD  condition,  six  from  LD12:12 
(four  mice  were  recorded  in  both  experimental  conditions) .   The 
data  within  each  section  were  then  "folded"  to  align  the  cycles 
along  circadian  phase  (LD12:12  sections  were  aligned  with  clock 
time,  DD  sections  were  aligned  using  the  onset  of  wheel-running 
activity  as  the  reference  phase) .   Data  were  averaged  across 
cycles  and  the  average  was  plotted  as  an  estimate  of  probability 
of  total  wake  or  total  sleep  as  a  function  of  circadian  phase  (7). 
Local  serial  correlation  (normalized  joint  second-order  prob- 
ability) was  estimated  across  cycles  at  each  circadian  phase  (8) . 
This  treatment  was  repeated  at  all  lags  in  the  range  0  -  500 
minutes  at  steps  of  one  minute. 

The  graphs  in  Figure  1  depict  mean  LSC  (averaged  across  all 
phases  of  the  circadian  cycle)  for  each  20-30  day  sample  plotted 
as  a  function  of  lag.   Note  that  LSC,  averaged  over  the  six 
LD12:12  and  eight  DD  sections,  is  significantly  positive  (i.e. 
greater  than  0.10)  for  all  lags  out  to  approximately  2  3  minutes. 
Most  important,  it  is  also  apparent  from  Figure  1  that  mean  LSC 
is  a  monotonically  decreasing  function  of  lag;  there  are  no  re- 
peated peaks  and  troughs  either  in  the  mean  correlation  curves 
or  in  any  individual  correlation  curve  that  might  reflect  the 
contribution  of  a  rhythmic  process  having  a  stable  period  in 
this  range. 


For  purposes  of  simplification,  only  the  representative 
correlational  results  obtained  at  a  lag  of  12  minutes  were  used 
to  illustrate  the  subsequent  analysis  of  phase-dependence  of 
LSC  (Fig.  2).   Figure  2  presents  the  LSC  functions  at  12  minute 
lag  for  wakefulness  from  a  representative  DD  section.   The 
probability  function  (marginal  probability  of  wakefulness)  is 
plotted  below  the  correlation  function.   Note  that  the  probability 
of  sleep  should  be  equal  to  one  minus  the  probability  of  wake 
and  the  correlation  functions  of  the  individual  halves  of  such 
a  binary  process  should  be  identical  (9).   In  actuality,  small 
imperfections  in  the  state  scoring  process  (occasional  unscored 
or  missing  data)  produce  some  deviation  from  the  expected  results. 
However,  these  differences  are  sufficiently  insignificant  to 
allow  presentation  of  the  data  for  wake  as  representative  of  the 
entire  process. 

The  figure  demonstrates  that,  although  mean  LSC  is  signifi- 
cantly positive  at  the  12  minute  lag,  the  value  of  the  LSC 
function  is  not  consistent  across  the  entire  circadian  cycle. 
Instead,  serial  correlation  is  maximal  at  phases  of  the  circadian 
cycle  corresponding  to  maximal  probability  of  wakefulness.   This 
specific  phase  dependence  of  the  local  correlation  structure  is 
a  consistent  feature  for  all  mice  observed,  though  the  relationship 
is  stronger  for  sections  recorded  in  DD  than  for  those  in  LD12:12 
(Table  1),  (10).   Additionally,  we  found  the  relationship  of 
correlation  and  phase  to  be  independent  of  lag  for  a  wide  range 
of  lags  in  the  ultradian  range  (10  seconds  to  20  minutes).   (As 
mentioned,  mean  LSC  approaches  zero  beyond  lags  of  20-25  minutes 
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TABLE  1 


Local  serial  correlation  (LSC)  of  wakefulness  at  circadian  maxima 
and  minima.   Each  entry  represents  the  average  of  the  90  LSC 
values  in  the  90  degree  phase  windows  encompassing  the  maximum 
and  minimum  probabilities  of  wakefulness.   A  t-test  was  used  to 
evaluate  the  difference  between  mean  correlation  in  the  two  phase 
windows . 


Mouse  and 
Condition 


MAX    P(WAKE)    MIN    P(WAKE) 


ALSC 


t-test 


PHASE 


LSC 


PHASE 


LSC 


LD01 
LD02 
LD0  3 
LD04 
LD05 
LD06 


DD01 
DD02 
DD03 
DD04 
DD05 
DD06 
DD07 
DD08 


14° 
2° 

11° 
5° 
5° 
9° 


7° 
6" 
2° 
7° 
2" 

16° 
2° 

13° 


656 
240 
417 
434 
298 
338 


626 
242 
453 
398 
355 
418 
299 
176 


264° 
256° 
263° 
212° 
226° 
245° 


241° 
266° 
252° 
270° 
252° 
270° 
270° 
260° 


.160 

.496 

.164 

.076 

X/S_=3. 556 

.189 

.228 

X 

.180 

.224 

.254 
.074 

(p    *     .02) 

.121 

.217 

2> 

1.345 

x= 

.224 

sx= 

.155 

s  = 

.063 

X 

.287 

.339 

.111 

.131 

X/S    =5.302 

.060 

.  393 

X 

.253 
.124 

.145 
.231 

(p    >     .01) 

.071 

.347 

.162 

.137 

.078 

.098 

£x=1.821 
X=  .228 


X 


115 


S_=  .043 
X 


and  significant  correlation  no  longer  exists  at  any  phase  of  the 
circadian  cycle.) 

The  consistency  of  positive  correlation  across  lags  and  the 
absence  of  significant  periodic  peaks  and  troughs  in  mean  LSC 
argue  against  a  stable  ultradian  rhythm  as  the  principal  source 
of  local  correlation.   Further,  at  no  lag  in  this  10  second  to 
20  minute  range  was  local  serial  correlation  found  to  be  sig- 
nificantly negative  for  any  circadian  phase.   Thus  it  is  unlikely 
that  the  observed  circadian  changes  in  correlation  are  the  result 
of  either  imprecise  estimation  of  the  ultradian  period  or  an 
ultradian  oscillator  whose  period  changes  stably  with  circadian 
phase.   Although  these  findings  alone  do  not  preclude  correlation 
arising  from  an  ultradian  oscillator,  they  do  suggest  that  such 
a  contribution  is  minor  relative  to  the  significant,  general 
local  serial  correlation  that  was  observed.   To  account  for 
all  observed  features  of  the  correlation  structure,  the  variability 
of  the  ultradian  "oscillator"  would  be  so  large  as  to  defy  the 
definition  of  the  term. 

A  more  plausible  explanation  for  the  observed  correlation 
structure  is  that  the  phase  of  the  circadian  cycle  corresponding 
to  peak  wakefulness  is  associated  with  a  simultaneous  increase 
in  mean  bout  length  for  both  states  and/or  an  increase  in  the 
serial  correlation  of  the  lengths  of  successive  bouts  of  a  given 
state.   That  is,  the  mechanisms  controlling  the  expression  of 
sleep  and  wakefulness  might  be  altered  at  this  phase  of  the 
circadian  cycle  such  that  the  lengths  of  successive  bouts  of  a 
given  state  are  correlated  to  a  greater  degree  than  at  other 


10 


circadian  phases  (11) .   Such  a  model  could  account  for  both  the 
dependence  of  serial  correlation  on  circadian  phase,  and  the 
occurrence  of  significant  serial  correlation  at  all  lags  in 
this  range.   Further,  serial  correlation  arising  from  such  a 
source  can  result  in  the  appearance  of  rhythmicity  in  frequency 
spectra  or  autocorrelations  (12),  a  phenomenon  that  may  account 
for  a  number  of  reports  of  ultradian  rhythmicity  in  physiology. 

Lastly,  it  is  clear  that  the  absence  of  evidence  for  ultradian 
rhythmicity  in  sleep  and  wakefulness  does  not  completely  preclude 
such  periodicity  in  the  REM  and  NREM  states.   However  preliminary 
results  using  the  same  analysis  techniques  presented  here  show  no 
evidence  of  ultradian  periodicity  in  REM  and  NREM  states  in  the 
mouse  (13).   To  remain  viable,  therefore,  it  would  now  seem 
that  models  proposing  ultradian  periodicity  must  either  abandon 
the  extension  to  a  rest-activity  cycle  as  proposed  by  Kleitman 
or  invoke  a  phylogenetic  origin  for  the  rhythm  after  the  evolution 
of  the  mouse.   Given  other  similarities  of  mammalian  sleep-wake 
organization,  this  latter  assumption  is  difficult. 
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The  local  serial  correlation  is  estimated  as 


f*^  ^^'^j^^V^S^f^^^V^S'^S^ 


1/2 


where 


111-  _L 

5(V*£+j>  =  Jn  X(k,£)X(k,£+j)/ 


m 


is  an  estimate  of  the  joint  probability  function. 


For  a  binary  time  series,  if   P    and   X.,  0.  =  1   correspon 

/\ 

to  state  1   and   P~   and   X ,,  p,  =  0   correspond  to  state  0, 


d 


then  we  can  compute   P..  ($  )  ,  P  (<j)  .,<})..  )   and   r,  (  £ , j  )   a 

7  and  8  previously. 
It  follows  that: 


s    in 


in- j. 


b)       VWj1    "    „L     ^-X(k,i)Hl-X(k(!l+j))/ 


111 


=        1       " 


m-1 

£     x(k  n  /  m! 

k=0        (     '     '  } 


m-1 

I     x 

k=0 


(k,£+j)   / 


rn 


m-1 

klQ    X(k,£)X(k,£+j)   /  mj 
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c)       rQ(£,j) 
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Iv^d-VV)  5o(*£+j>f1-5o(*£+j>5]1 


/2 


(1   -  pi(4>£)    -  Pl(^+j)    +  p^VW 
-    (i-p1((Dji))(i-p1((J)        )) 


[(l-P1(^))(P1(^)[l-P1(^+j))P1(*£+j) 
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P1(  +  &'W     -    WVW 


P^^jfl-P^^))     P1(^+j)(l-P1(^+.))_ 


1/2 


r1(^,j) 


10.   The  apparent  weakening  of  phase-dependence  of  correlation 
observed  in  LD12:12  sections  is  not  thought  to  represent  a 
true  effect  of  lighting  condition.   Rather  it  is  believed  to 
be  the  result  of  reduced  variation  in  the  phase  alignment  of 
these  sections  resulting  in  a  larger  portion  of  the  cycle 
in  which  the  marginal  probability  of  wake  is  near  or  equal 
to  one  (p (sleep)  =  O) .   The  correlation  behaves  erratically 
when   P   is  very  near  to  one  and  this  effect  results  in  a 
reduced  average  correlation  across  the  90  degree  window 
encompassing  maximum  probability  of  wake.   For  DD  sections, 
the  increased  phase  variability  of  the  free  running  rhythm 
results  in  fewer  points  where  the  marginal  probability  is  very 
close  to  one  and  consequently  the  correlation  function  is 
better  behaved. 
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